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Abstract: We use holographic techniques to study meson quasiparticles moving through a 
thermal plasma in N = 2 super- Yang-Mills theory, with gauge group SU(N C ) and coupled to 
N { flavours of fundamental matter. This holographic approach reliably describes the system 
at large N c , large 't Hooft coupling and N f /N c <C 1. The meson states are destabilized by 
introducing a small quark density n q . Spectral functions are used to examine the dispersion 
relations of these quasiparticles. In a low-momentum regime, the quasiparticles approach a 
limiting velocity which can be significantly less than the speed of light. In this regime, the 
widths of the quasiparticles also rise dramatically as their momentum approaches a critical 
value q crit . While the spectral functions do not display isolated resonances for q > q crit , the 
dispersion relations can be extended into this high-momentum regime by studying the dual 
quasinormal modes. A preliminary qualitative analysis of these modes suggests that the group 
velocity rises to the speed of light for q S> q crit . 
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1. Introduction 



A large class of strongly coupled gauge theories can be studied using the gauge/gravity duality 
[1, 2]. The gauge theories that are currently amenable to study with holographic techniques 
are very different from real world QCD, e.g., current calculations are restricted to large N c and 
large 't Hooft coupling. However, this approach has still proven to be a fruitful framework 
with which to gain new insights into the strongly coupled quark-gluon plasma - see, e.g., 
[3,4]. 

With this aim in mind, holographic techniques have been applied to study the thermal 
properties of N { flavours of fundamental matter in N = 2 SU(N C ) super- Yang-Mills (SYM) 
in a quenched approximation (i.e., N { <C N c ) [5, 6, 7, 8]. The gravity dual for this field 
theory consists of N { probe D7-branes in the black hole background generated by N c D3- 
branes. In this system, the fundamental matter generically undergoes a first order phase 
transition at some temperature T fun . The low-temperature phase of the theory is described 
by 'Minkowski embeddings' of the probe branes (see fig. 1) in which the branes sit entirely 
outside the black hole [6, 7]. In this phase, the meson spectrum is discrete and exhibits a 
mass gap. Above the critical temperature T fun , the branes are characterised by 'black hole' 
embeddings which extend through the event horizon. In this phase, the meson spectrum is 
continuous and gapless [7, 9, 10]. Thus, this large- N c , strong coupling phase transition is 
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associated with the dissociation of the mesons. It has been suggested that this physics is in 
qualitative agreement with that of heavy quarkonium in QCD [7]. Studies from lattice QCD 
[11] suggest that such mesons survive the deconfinement phase transition at T dec ~ 175 MeV 
and remain as relatively well-defined resonances up to temperatures of 2 — 3T dec . 



Minkowski embedding Black hole embedding 




Figure 1: Possible embeddings for probe D7-branes in the background black hole geometry of D3- 
branes. 

Another interesting feature that was discovered for the holographic mesons in the low 
temperature phase is that moving through the thermal plasma, at large momentum they 
approach a limiting velocity with <y lim < 1 [7, 12]. On the gravity side, this can be understood 
as the dual excitations travelling at the local speed of light near the minimum radius reached 
by the D7-branes. Hence in the field theory, one finds the velocity 

dx 



dt 



Vlim = -77 = V~9tt/gxx (1-1) 



P=Pn 



Because of the redshift near the black hole horizon, this yields a result which can be much 
less than one. As discussed in [12], this result (1.1) can be re-expressed as 



4 \M q/ 

We must note that the limiting speed does not actually represent a 'speed limit' for the mesons. 
That is, a careful analysis show that the group velocity actually approaches t> lim from above 
and so the maximum group velocity is actually slightly larger than this asymptotic value [12]. 
In any event, the effect that t> lim < 1 is somewhat surprising as one's naive intuition would be 
that a meson traveling through the thermal plasma would eventually reach the speed of light 
if the energy /momentum is increased to arbitrarily large values. Hence it appears that this 
new limiting velocity is a consequence of strong coupling. 
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As mentioned above in the low temperature phase, the holographic mesons are stable or 
rather their widths are suppressed by 1/N C . Of course, this stands in contrast with the heavy- 
quark mesons studied in thermal lattice QCD [11] which are quite broad states. In the black 
hole phase of the holographic model, the meson excitations are readily absorbed by the black 
hole horizon. The quasinormal frequencies of these excitations typically have Im(w) ~ Re(cj) 
[10] and so the corresponding spectral functions do not reveal any quasiparticles [9]. However, 
black hole embeddings can also arise at low temperatures when a nonzero quark density n q 
is introduced 1 - in fact, these are the only physically consistent embeddings in this situation 
[13]. In these backgrounds, the meson states are again unstable but the width of these states 
can be tuned by varying n q [14, 15]. Essentially, with a small ra q , the D7-brane extends down 
to the horizon with a narrow neck (as illustrated in fig. 2) and so the absorption of the meson 
excitations is limited by the small effective horizon area of the worldvolume metric. Hence the 
corresponding spectral functions display clear quasiparticle states with narrow widths [15]. 
In this paper, we begin a study of the dispersion relations of these quasiparticle states. 




Figure 2: D-brane embedding with radial electric field is necessarily a black hole embedding but the 
width of the throat is tuneable. 

An overview of the paper is as follows: in section 2 we introduce the D3/D7-brane 
framework and review the D7-brane embeddings and thermodynamics. In section 3, we 
turn to the calculation of the spectral function for various meson operators. By following 
the position and shape of the quasiparticle resonances in the spectral function with growing 
momentum, q, we estimate the dispersion relations for the low-lying resonances. With this 
approach, the quasiparticles are found to approach the same limiting velocity found for the 
case of stable mesons [7, 12] 2 . Another interesting phenomena is that the widths show a 
dramatic increase as the momentum approaches a critical value q crit . Section 4 presents a 

1 Since the the underlying theory here is supersymmetric, the quark density we consider arises from a 
hypermultiplet with both fermions and bosons. 
2 See also [16]. 
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qualitative discussion of the quasinormal frequencies using our intuition derived from casting 
the relevant radial equation in the form of an effective Schrodinger equation. This framework 
gives an alternative point of view from which to understand the effects noted above. Above 
q crit , the spectral functions do not exhibit isolated quasiparticle resonances and so in section 
5, we use the Schrodinger framework to consider the behaviour of the quasinormal modes in 
the high-momentum regime q > q crit . In previous sections, we defined the limiting velocity 
Vnm by observing the (real part of the) dispersion relation approached a straight line for large 
q but in the regime q < q crit . Here our qualitative analysis suggests that this definition fails 
in the regime q » q crit . Rather the behaviour of the dispersion relation changes such that 
ultimately it approaches an asymptotic slope of one in this very high momentum regime. 
In section 6, we discuss our results and make a few observations about future directions. 
Appendix A provides some details about a WKB calculation of the quasinormal frequencies 
using the Schrodinger framework considered in section 4. 

2. Holographic Framework 

Following [6, 7], we write the background metric for N c black D3-branes in the decoupling 
limit as 



f 2 

-^dt 2 + fdx 2 

f 



+ — [dp 2 + P 2 dn 2 5 ] , (2.i) 



(2.3) 



2 V L J 

where p is a dimensionless coordinate and 

f(p) = l-^i, f{p) = l + j A , L 4 = 4Trg s Njl (2.2) 
This metric possesses a horizon at p = 1 with temperature 

T = -^L 

ttL 2 

which is identified with the temperature of the dual J\f = 4 SYM theory. Further the coor- 
dinates {t, x} are identified with the coordinates of the gauge theory. The string coupling 
constant is related to the SYM 't Hooft coupling constant through 3 

A = g 2 M N c = 2ng s N c . (2.4) 

The background D3-brane solution also has a (constant) dilaton and a Ramond-Ramond field, 
whose precise form are not needed in the following. 

Introducing N t D7-branes into the geometry above corresponds to coupling N { fundamen- 
tal hypermultiplets to the original SYM theory [17]. Before the decoupling limit, the branes 
are oriented in the following array: 

0123456789 
D3: x x x x (2.5) 
D7: xxxxxxxx 



3 Note that we are using the standard D-brane convention here which differs from that of the usual quantum 
field theory literature. As explained in appendix D of [7], A = Aqft/2. 
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This configuration is supersymmetric at zero temperature, which ensures stability of the 
system. After the decoupling limit, the D7 branes wrap an S 3 inside the S 5 of the background 
geometry. Adapting the S 5 coordinates to this embedding, we write 

= d6 2 + sin 2 6 dQ 2 3 + cos 2 6 # 2 (2.6) 

in (2.1). Defining x = cos#, we describe the D7-brane embedding as: <j> = 0, x = x(p)- 

The derivation of the equations of motion for the D7-brane profile x(p) an d the gauge 
field on their worldvolume At was discussed in [13]. Here we will review a few salient points 
and refer the interested reader to [13] for more details. The DBI action of the D7-branes may 
be written 



J D7 = —N f T D7 J dtd 3 xdpdn 3 ^^ff(l - X 2 )^l- X 2 + P 2 X 2 _ _ X 2)i t 2 ? (2 . 7) 

where the dot denotes differentiation with respect to p and we have introduced the dimen- 
sionless gauge field A t [13]. The asymptotic form is determined by the gauge field's equation 
of motion (eqn. (2.11) in [13]) as 

~ 2n£ 2 d , , 

A t = i Mq -_ + ..., 2.8 

u p 2 

where the constant /U q is the quark chemical potential. The dimensionless constant d is related 
to the vacuum expectation value of the quark number density operator with 

n, = ^N t N c V\T 3 d. (2.9) 

The equation of motion for % (eqn. (2.17) in [13]) implies that the D7-brane profile behaves 
asymptotically as 

m c . 
X = - + - + •••, (2.10) 
P P 

where the dimensionless constants m and c are proportional to the quark mass and condensate, 
respectively [6, 7]. In particular, m = M/T where 

2M Q M EaD 

M = ^- = ^^- (2.11) 

V\ 2vr v ; 

is (up to a factor) the meson mass gap M gap at zero temperature [18]. 

As only d p A t enters the action (2.7), J is a conserved integral of motion and gauge field 
equation (Gauss' law) yields 



3 p A t = 2d - X 2 + p2 * 2 . (2.12) 

f{i-x 2 )[p & P(i-x 2 ? + m 
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Substituting this expression into the equation of motion for the profile x then gives 



a i P 5 ff(i-x 2 )x 

where A is given by 



2.1 Black hole embeddings 



//>/A =lmi-X 2 + P 2 X 2 )- P 2 X 2 ] ■ (2.13) 



VA y /l- X 2+p2 x 2 



A = g. 1 {L — X > _ . (2.14) 

P 6 / 3 (l-X 2 ) 3 + 8^ 2 



A key point for the following analysis is that if n q ^ (i.e., d ^ 0) then the only physically 
consistent embeddings for the D7-branes are black hole embeddings [13, 14]. Simply stated, 
a nonzero density of quarks is dual to a worldvolume electric field, i.e., a nonvanishing A%. In 
turn, this electric field can be interpreted as a finite number of (fundamental) strings dissolved 
in the probe D7-branes. Since these strings cannot simply terminate, it is not possible for 
the D7-branes to close off smoothly above the horizon. 

Of course, the defining feature of the black hole embeddings is that the probe D7-branes 
reach the event horizon at p = 1. We can only solve the profile's equation of motion (2.13) 
using numerical techniques. Generally we integrate out from the horizon with the following 
boundary conditions 

x ( P = l) = Xo , d pX (p=l) = 0. (2.15) 
The series expansion of x around p = 1 takes the form 4 

X = Xo ~ J^ (1 ~ X °C 3 > - !) 2 ( 2 ~P) + °«P ~ !) 4 ) • ( 2 - 16 ) 
4(d 2 + (1 - xlf) 

Thus as xo approaches 1 (from below) with d fixed, the profile develops a long narrow throat 
extending out from p = 1 - as illustrated in fig. 3. 

The black hole embeddings were so named because the induced geometry on the D7-brane 
worldvolume is a black hole geometry. The area of the induced horizon, which is proportional 
to (1 — Xo) 3//2 ' controls the lifetime of excitations on the brane [9]. As discussed below, tuning 
Xo close to 1 will allow us to produce long-lived quasiparticles. However, in the following, 
we will want to compare the lifetimes of various quasiparticles while keeping the quark mass 
fixed. Recall that the latter is determined from the asymptotic form of the profile (2.10) at 
large p. Implicitly then with our numerical approach, the constants m and c in (2.10) are 
functions of both xo and d. Hence adding the parameter d is the key to allowing us to vary 
Xo (i.e., the quasiparticle lifetimes) while holding m (i.e., the quark mass) fixed. 

We note that the zero-temperature limit of these configurations is a subtle one [19]. 



4 In the numerics, we used this expansion to specify the boundary conditions slightly away from p = f . 
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J 1 1 1 1 1 1 1 1 1 1 1 1 1 ,_ p 

2 4 6 8 



Figure 3: The profile x as a function of p for different values of xo, d with fixed m = 4.56. To keep 
m fixed, d is decreased as xo increases . The black line represents the profile for the corresponding 
Minkowski embedding with d = 0. 



2.2 Minkowski embeddings 



We comment briefly on Minkowski embeddings with d = with reference to the limiting 
velocity discussed in the introduction. To describe these embeddings, we make the following 
coordinate transformation [7] 



p 2 = r 2 + R 2 , X = R/P, 
whereby the equation of motion for the profile R{r) becomes 



(2.17) 



r 3 ffd r R 



8r s R 



(r 2 + R 2 )' 



■y/l + (d r R) 2 



(2.18) 



where / and / are defined as in (2.2) but now expressed in terms of r and R using (2.17). 
The boundary conditions specified at r = (i.e., \ = 1) are: R = Ro,d r R = with Rq > 1. 
Near the axis r = 0, the embedding has an expansion. 



R — Rn + 



Rq(Rq 



I) 



+ 



5Rl 6 + 5R% 



3RURI 



1) 



0(r e 



(2.19) 



which, as expected, clearly illustrates that r = is the point of closest approach between the 
D7-brane and the horizon. As commented above, the meson excitations on the Minkowski 
embeddings were found to approach a limiting velocity at large momentum [7]. This velocity 
is given by 

= / = 1 - l/fi 4 

p=p mm f r=0 1 + 1 / R 



9tt 

Qxx 



4 - 



(2.20) 
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3. Meson spectral functions 



A fruitful application of gauge/gravity duality techniques has proven to be the study of 
thermal properties of the strongly coupled gauge theories through their finite-temperature 
correlation functions [20]. In a holographic framework, the spectral functions are typically 
easier to compute than the full correlators and the poles (and associated residues) of the corre- 
lators are still reflected in the corresponding spectral functions. According to the holographic 
dictionary, the poles are determined by the quasinormal spectrum of a dual bulk field fluc- 
tuations, whose study is a technically challenging problem [21] — we begin to address these 
calculations in sections 4 and 5. In contrast, the spectral function is given by the imaginary 
part of the retarded correlator, 

*(w,q) = -2LnG%;,q), (3.1) 

and is determined by appealing to standard calculations of bulk field correlators [22]. 

The retarded correlators have poles in the lower half-plane of complex frequency, which 
we may assume have the form 

G R ~ , (3.2) 

u - (l(q, a) + iT(q, a) v ' 

where a represents any additional relevant parameters, e.g., the temperature or quark density. 
From such a pole, the spectral function receives a contribution 

2 AT 

""'- (■-nr + r- - (3 ' 3) 

Thus in the vicinity of uo = 0,, the spectral function has a peak characterized by a width V. 
A quasiparticle interpretation can be given to the peak if it satisfies the Landau criterion: 
r <C ft. Hence the masses and lifetimes of quasiparticles can be extracted from the holographic 
spectral functions. 

The spectral function H(u;) also has a characteristic form in the 'high-frequency' limit. 
This behaviour is determined by the leading short-distance singularity 



A 

(t2"x^0^ V "' ~ ~" \t 2 - X : 



lim J0(t,x)0(O)) = - 9 „ 2 . A +---, (3.4) 



where A denotes the dimension of the operator O and A is a dimensionless constant. A 
Fourier transform then leads to the following contribution to the spectral function 



~ A{lo 2 - q 2 )*- 2 . (3.5) 

Note that this high-frequency tail is Lorentz invariant and shows no indications of a lim- 
iting velocity. This should be expected as it describes the very high-energy/short-distance 
behaviour, which is independent of temperature [9]. 
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In the present study on the gravity side, the D7-brane embeddings extend through the 
event horizon of the AdSs black hole, which describes the theory at finite temperature. As 
indicated by the name 'black hole embedding', the metric induced on the worldvolume of the 
D7-branes is itself a black hole. Even though the latter geometry does not obey Einstein's 
equations, the analysis of the hydrodynamic physics found previously for bulk fields, e.g., [20], 
is readily transferred to the worldvolume fields on the D7-brane [9] and hence we can examine 
the spectral function for various mesonic operators following the techniques introduced in [23]. 

Here, the area of the horizon induced in the worldvolume geometry controls how quickly 
excitations on the D7-branes are absorbed by the black hole. This absorption rate then de- 
termines the lifetime or width of the corresponding quasiparticles in the dual gauge theory. 
For a black hole embedding with a large horizon area {e.g., the typical situation in the high 
temperature phase when d = 0), one expects to find that r ~ 0, [10]. Further the corre- 
sponding spectral functions are essentially featureless beyond exhibiting the characteristic 
high-frequency tail (3.5). However, following the discussion in the previous section, by tuning 
Xo close to 1, the horizon area and consequently the quasiparticle widths shrink. Hence this 
tuning can bring us into a regime where the spectral functions display distinct peaks with 
r <C ^ [15]. We re- iterate that having the freedom to independently tune both d and xo is 
the key to producing these small widths while holding m (i.e., the quark mass) fixed. 

In the following, we calculate the spectral densities for the flavour current which is 
dual to the worldvolume gauge field on the supergravity side [13]. These calculations are 
a simple extension to finite spatial momentum of the spectral function calculations appearing 
in [15]. We begin by writing the full gauge field as 5 



A^ denotes a fluctuation. For simplicity, we assume here that A^ only depends on p and 
the Minkowski coordinates (and not the internal coordinates on the S 13 ). To determine the 
linearized equations of motion for the fluctuations, we expand the DBI action to quadratic 
order. The resulting gauge field lagrangian is 



where Q = g + F being the sum of the background metric and background electromagnetic 
field. Note that the last term does not vanish since Q is no longer a symmetric matrix in 
general - see comments below, however. The linearized equations of motion for A are then 



4(p,x) = ^I t (p)+^(p,x). (3.6) 
Here A t denotes the background gauge field (2.12) producing the finite quark density, while 




(3.7) 




(3.8) 



Next we consider the Fourier transform of A 




(3.9) 



implicitly we are working with the dimensionless gauge field scaled as in [13]. 
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If we restrict the momentum vector to (uj,q, 0,0), the components of the electric field are 
given by 

E x = ujA x + qA , E V)Z = ojA VjZ . (3.10) 

In the analysis which follows in this section, we will focus on the transverse fields E VjZ 
both of which we will denote by Et- The analysis for the longitudinal component E x is more 
involved on two counts. First, it is only for this mode that the second term in the lagrangian 
(3.7) contributes to the equation of motion. The second complication is that, in fact, the 
quadratic lagrangian above is inadequate to describe these modes. This is because with the 
background field F pt , these longitudinal modes mix with the scalar fluctuations 59. Hence we 
leave the complete analysis of these coupled modes for the future. 6 

The linearized equation of motion for Et can be written as 



^2 



with 



F = ^L - (1 -* 2)2 , G = ^V 1-X 2 + PV 2 (1-X 2 )- (3.12) 

VA y/l - X 2 + P 2 X' 2 P f 

Here we have also introduced to = uj/(2ttT) and q = q/(2nT). Further recall that A is defined 
in (2.14) and goes to unity when d = 0. The retarded Green function for Et is given by 

F d p E T 



C R NfN c T 2 



lo 2 Et 



(3.13) 



Now the spectral function is defined as 

k) = -2 to 2 ImG R , (3.14) 

where the extra factor of J 2 is introduced to produce the spectral function for A y>z - compare 
to (3.1). Defining 

F d E T 

s = rf> ( 315 » 

it follows from (3.11) that 



_2 / f2 J2 



W + T s + G {Kp-i?) =0 - (3 - 16) 

Thus we have 

*(u,q) = I^«> 2 I m 5(p^oo). (3.17) 

We solved the above equation numerically taking Et(p) = {p — l)~ tto e(p) where e(p) is 
regular at the horizon with the boundary conditions: e(l) = 1, d p e(l) = ivo/2. This implies 

5(1) = -4i^(l- X §) 3 + dV«>. 

6 These complications are not evident from the analysis in [15] due to simplifications at q = 0, for e.g., 86 
only mixes with Aq. 
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Figure 4: Spectral functions for m — 4.56 and q = 0. The vertical black lines correspond the 6- 
functions appearing for the Minkowski embedding. The latter appear above at v» = 4.56, 7.89, 11.16, 
respectively. The green line is for d = 0.06, \o = 0.99755, red for d = 0.15, \o = 0.99394 and blue for 
d= 0.25, xo = 0.99. 



Fig. 4 illustrates the behaviour of the spectral functions at q = - this is reproducing the 
results given in [15]. As discussed above, to construct the background D7-brane embedding, 
we are varying the variables xo an d d while keeping fixed the asymptotic mass (or m). For 
the cases illustrated in fig. 4, we fixed m = 4.56 (to within about 0.25% accuracy). These 
examples explicitly show that as the quark density or induced horizon area increases, the 
widths of the quasiparticles increases while their positions remains essentially fixed. The 
figure also shows the positions of the corresponding mesons on the Minkowski embedding 
with m = 4.56 (and d = 0) - these masses were determined numerically, as described in [7]. 

3.1 Dispersion Relations 

Now to study the dispersion relations, we consider the spectral functions with finite spatial 
momentum q. We focus our attention on the first few peaks as these lowest lying resonances 
are the most prominent at q = 0. At least for a certain range of q, the spectral functions allow 
us to reconstruct both f2(q) and T(q) for these quasiparticles. As might be expected, we will 
find that beyond a certain q crit , the corresponding poles (3.2) have moved too far away from 

In all our plots in the following, the spectral function 'K is in units of f / — . 
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the real axis to allow us to identify individual resonances in the spectral functions. We will 
also be able to estimate the residue A(q) in the regime q < q crit . 

Fig. 5 illustrates the typical behaviour of a spectral function as the momentum is in- 
creased. In general, one can observe that the quasiparticle peaks are moving to larger values 
of to as q increases. Similarly, one would say that increasing q causes the individual peaks 
to become diminished while the background increases. In particular, the high frequency tail 
(3.5) in the present case is 

*(w) = 4vr (n> 2 - q 2 ) (3.18) 

with A = 3 [9]. In the fig. 5, we see that fairly quickly as the momentum is increased, the 
spectral function becomes well modeled by this tail alone. In fact in the last two plots with 
q = 100 and 140, the rise in the spectral function cannot be distinguished from this tail at 
the scale used there. 
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Figure 5: Spectral function versus no with d = 0.25, xo = 0.99, m = 4.56 and increasing q, as 
indicated. There are no visible peaks above q = 100. In each plot, the dashed grey line indicates the 
high frequency tail (3.18) for vo > q. 



Energies, Q(q): To make our discussion more quantitative, we followed the positions fij(q) 
of the (first few) peaks as q increases for three different limiting velocities. The three cases 
are summarized in the following table: 







m 


d 


Xo 


n 


r 


I 


0.995 


4.56 


0.25 


0.99 


4.52 


0.072 


II 


0.651 


1.50 


0.0005 


0.9995 


1.37 


0.002 


III 


0.343 


1.32 


0.0001 


0.99975 


1.03 


0.001 



(3.19) 
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10 20 30 40 50 2 4 6 8 10 12 14 

Figure 6: The plot on the left shows the positions of the first three peaks for case I (see table (3.19)) 
as a function of q. Within our accuracy, these curves all asymptote to roughly to = 0.995q, which 
corresponds to the black line. The plot on the right displays similar dispersion curves for all three 
cases. The blue, red and green lines correspond to cases I, II and III, respectively. 

The columns labelled f^o and Tq in the table above indicate the position and width of the 
first peak at q = 0, respectively, for each case. In each of these cases, we chose Tq/Qq <C 1 as 
this allowed us to follow the first few peaks to much higher momenta. As we will see below, 
if one begins with broader peaks at q = then they dissipate more quickly as q is increased. 
As described above at (3.3), the positions of such the peaks 8 Slj(q) should correspond to the 
real part of the position of a pole in the corresponding retarded correlator. 9 

The resulting plots of f2i(q) for these three examples are displayed in fig. 6. Note that 
fig. 6 contains two plots which differ primarily in the horizontal scale for q because while the 
peaks could be followed out to q « 100 for case I (u lim = .995), they disappeared around 
q 10 in the last two cases - see discussion below. In all three cases, the curves fij(q) appear 
asymptote to the straight-line form: 

fii(q) = u lim q + a; + 0(l/q) (3.20) 

found analytically for the Minkowski embeddings in [12]. Further, as is implicit in our no- 
tation, the asymptotic value of the slope <9qf2j matches very well the limiting velocity v lim 
calculated from the corresponding Minkowski embedding with the same value of m. For 
^um = 0.995, this is best illustrated in the left panel of fig. 6. Fig. 7 shows the first peak for 
the three different cases fit with a straight asymptote with slope v Um , as given in (3.19). This 
figure makes clear that the two cases with slower limiting velocities approach their asymptotic 
behaviour (3.20) more quickly. We should comment that while the limiting velocities in (3.19) 
gave a good fit to the asymptotic behaviour of these curves, the accuracy of our numerical 
calculations was limited. The numerical error is estimated as follows: Using t> lim ~ J ) , we 

s In the following, the subscript is i = 1,2, 3 to denote the first three resonances in the spectral function. 
Further we identify the position of the resonance with the condition d^iTi = here. 

9 In the following discussion, we denote the position of the poles in the complex to-plane as Q, — iT. Hence 
this usage differs from that in (3.2) or (3.3) by a factor of 2nT. 
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q 



Figure 7: The points above correspond to positions of the first peak from the spectral function 
for each of the cases in table (3.19). The corresponding straight lines, which match the asymptotic 
behaviour in each case, have slopes v lim = 0.995, 0.651 and 0.343, respectively. The dashed black line 
has slope one, corresponding to null four-momenta: f2 2 — q 2 = 0. 



have Av llm = 8AR -S-v liin . Using AR = ±0.01, this gives Av lim « ±0.00004, ±0.008, ±0.01 
for v Vlm = 0.995,0.651,0.343 respectively. Hence we are being slightly extravagant in quoting 
three significant figures for the latter two cases. We also note that we found that the constants 
cij did not seem to match well with the analytic expressions given in [12] for the Minkowski 
embeddings. 

Widths, r(q): Next we wish to examine how the widths of the quasiparticles evolve with 
increasing momentum. For this purpose, we list four new trial cases in the table below. Note 
that the initial widths Tq are tuned to be roughly equal in cases IV, V and VI, which were 
also chosen to match the same values of v lim as appeared in cases I, II and III in table (3.19). 
We have also chosen case VII with v nm = 0.651 but tuned so that the ratio Tq/Qq is roughly 
the same as in case IV with v lim = 0.995. 







m 


d 


Xo 


n 


r 


IV 


0.995 


4.56 


0.15 


0.99394 


4.52 


0.044 


V 


0.651 


1.50 


0.01 


0.99 


1.37 


0.04 


VI 


0.343 


1.32 


0.0033 


0.991 


1.03 


0.042 


VII 


0.651 


1.50 


0.004 


0.996 


1.37 


0.015 



(3.21) 



As described at (3.3), near the location of an isolated peak in the spectral function, we 
expect that the spectral function can be approximated as 



2.4T 

(t»-o) 2 + r 2 ' 



(3.22) 
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Prom this form, we can derive the following expression for the width of the peak (i.e., the 
imaginary part in the position of the corresponding pole): 



r 



_ 2 A 



(3.23) 



where ' denotes differentiation with respect to to. In the above expression, we define to = 0, 
as the point where djx,l& = 0. The results for the width of the first peak as given by this 
formula are shown in fig. 9 (for each of the four cases in table (3.21)). Note that each of these 
curves shows a dramatic increase in T as q increases from zero. This gives the first hint that 
we should find a certain maximum value of q beyond which the quasiparticles do not exist. 
Certainly this is observed, e.g., in fig. 5 where the peaks in the spectral function are simply 
washed out at large q. 

Note however, that if the spectral func- 
tion contains a significant background con- 
tribution near vo = O, then the above 
formula (3.23) will tend to overestimate 
the true value of T. In particular, we can 
' seek to improve the approximate form of 
the spectral function (3.22) in the vicin- 
ity of one of the peaks by also including 
the high frequency tail (3.18) 




Figure 8: A plot of the first peak in "H (red) and the 
corresponding c^'H (blue) for case V with q = 0. This 
illustrates the positions of the frequencies t»o, t»o,± and 
t» m> ± discussed in the text. <9t»^ nas be scaled by 
1/200 in this plot. 



2vir 



+ 47r(o> 2 - q 2 



(»-fi) 2 + r 2 

(3.24) 

Now we observe, however, that by tak- 
ing the third derivative, c^'X, the back- 
ground introduced by the tail in the above 
expression will be eliminated. Further 
the resulting expression will have zeros at i»o = ^ an d n>o,± = SI ± T and so in princi- 
ple, the zeros of can be used to determine T. However, rather than working directly 
with these zeros, we note that there are also a maximum and a minimum to either side of the 
central zero at £1. A straightforward calculation shows that these extrema occur at 



t» m ± = Q, ± 5 r 



with 5 



0.3249. 



(3.25) 



The above features are illustrated in fig. 8. We can use these expressions to estimate both the 
width of a given peak but also an error because our form (3.22) is not perfect. In particular, 
we have 



26 



AF = - K 



2m 



■o) 



(3.26) 
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Figure 9: log(r/rq =0 ) vs log(ci + 1). The continuous curves were obtained using the spectral plots 
and formula (3.23). The blue, red, green and brown lines correspond to cases IV, V, VI and VII in 
table (3.21), respectively. The points with 'error bars' were obtained using (3.26). The vertical dashed 
lines within the shaded regions indicate the location where the zero t» mi + in d^T^ is lifted. The left 
extreme of the shaded regions correspond to q when t»o+ is lifted while the right extreme corresponds 
to q when t»o is lifted. The solid vertical lines indicate where a point of inflexion occurs in the effective 
potential, as described in section 4. The arrows indicate where q = fi. 



The results for the width using these expressions are also shown in fig. 9. One can note that 
again there is dramatic rise in T as q increases, again hinting at a q crit . These results for T 
seem to agree fairly well with those originally derived with (3.23). However, this agreement 
is not as good for larger values of q where the quasiparticle peaks are washed out and AT is 
larger. In particular, as might be expected, (3.23) tends to give a larger width than (3.26) in 
this regime. 

In fact, the increase in the widths is only one of three effects leading to the dissolution 
of the quasiparticle peaks in the spectral functions. The second effect which we consider here 
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is the decreasing separation of neighbouring peaks or poles along the real axis. As shown in 
fig. 9, one finds that Ar from (3.26) also increases with increasing q. This is indicating a 
greater asymmetry in around the position of the peak xoq = CI. Of course, this is actually 
indicates that the form (3.24) is becoming a poor approximation for the spectral function in 
the vicinity of the peaks. Essentially the problem is that assuming that the peaks are isolated 
is incorrect in this regime. To illustrate this behaviour for case VI, we show the ratio of the 
separation of the first two peaks to the sum of their widths in fig. 10. Prom our plot, we see 
that this ratio approaches one indicating that this first peak is no longer isolated from the 
other quasiparticle peaks at higher rx>. The same behavior is observed for the other cases. 
For cases IV, V, VII we have for q = 0, the ratio log(ft 2 - ^i)/(f 2 + fi) = 2.54, 1.54,2.36 
respectively while the ratio near q where t» m ,+ becomes a point of inflection becomes .28(q = 
90), .40(q = 3.6), .36(q = 5) in the respective cases. 

The first peak which we are following then loses its shape due to the encroachment of 
the neighbouring peaks. The typical development as q increases is that first the zero at i»o,+ 
is lifted. Next the maximum at t» mj + becomes a point of inflection and then is subsequently 
lost. Finally, the zero at t»o itself lifted when it collides with t»o,- (and simultaneously with 
t» mj _). The region where the first peak is losing its shape in this way is indicated with the 
shading in fig. 9 for each of the four cases. This dissolution of the first peak is illustrated in 
detail for case V in fig. 11. As shown, first t»o,+ is lifted around q = 2.6. This is followed by 
the maximum t» mj + becoming a point of inflexion around q = 3.7 and finally t»o is itself lifted 
around q = 4.9. 

Beyond the point where t»o is lifted, there is no (simple) way to use the spectral function 
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Figure 10: Plot of log(f22 — f2i)/(r2 +Ti) for case VI. The central values arise from using T 1,2 = f 1,2 
while the upper and lower 'error' bars are determined with I\2 = Ti.2 ± AI\2- The ratio approaches 
one for larger q indicating that these resonances in the spectral function should not be considered to 
be isolated peaks. 
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Figure 11: The central plot shows the behaviour of T with q for case V as in fig. 9. For each point 
in this plot, we also show the corresponding spectral function X denoted by red and denoted by 
blue. This typical example illustrates how the peaks lose their form with increasing q. 

to follow the corresponding pole in the thermal correlator and certainly there is no sense in 
applying the concept of a quasiparticle. One might use the lifting of the zero t»o as defining 
the value q crit where the quasiparticles disappear. These momenta for each of the four cases 
in fig. 9 coincides to the extreme right of the corresponding shaded region. We return to these 
issues in the discussion in section 6. 

A clear trend that emerges from fig. 9 above is that the value of q crit increases with 
increasing u lim . This is already clear in comparing the results for cases IV, V and VI, all 
of which start with To — 0.04. However, while v Um is certainly one parameter that changes 
between these three cases, another parameter which distinguishes these three cases is Tq/^Iq. 
The latter increases significantly between these three cases while the limiting velocity is 
decreasing and one might imagine that if the quasiparticle peaks begin by being less well 
resolved, they should be washed out more quickly. Hence we also considered case VII for 
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Figure 12: The plot shows log(r2/Ti) for the first two peaks in case VI. For these two peaks, Tq 
takes the values 0.042,0.184 respectively. The first peak disappears at q ~ 3 while the second peak 
disappears at q« 2.72. 



which u Um = 0.651 and Tq/Qo — 0.01. So here the limiting velocity matches that in case V 
(for which To/^o — 0.03) but the ratio Tq/CIq was tuned to match roughly that in case IV 
(for which u lim = 0.995). In fig. 9, we see that the width for case VII grows with q in close 
accord with that for case V. We note though that our estimates of q crit for case VII were 
slightly larger than those of case V. In any event, these results seem to indicate that v Xim is 
the dominant factor in determining critical momentum where the quasiparticles disappear. 

For comparative purposes, we also show in fig. 12 the ratio of the widths of the first two 
peaks in case VI, as calculated with (3.26). As the plot shows, the width of the second peak 
increases more slowly than the first. For instance, r2(q)/Fi(q) ~ 4.3 at q = but only 1.3 at 
q = 2.2. This behaviour is in fact typical for all of the cases which we studied here. 

Residues, A(q): So far our discussion of the quasiparticle peaks has focussed on their energy 
O, and the width T, which correspond to the real and imaginary parts of the position of a pole 
in the thermal correlator (3.2). Here we briefly turn to the residue A, as this also determines 
the size of the peak corresponding to a given pole - e.g., for the canonical shape (3.3), the 
maximum value of spectral function is: ^(ro = f2) = 2A/T. Assuming the peaks have this 
canonical shape, the residue can also be determined in a number of different ways from the 
spectral functions but we focus on the third derivative of the spectral function following the 
above discussion for the form given in eq. (3.24). First, we can relate A to the slope of dm X 
at to = i»o : 

£ 5 

48 



Ai = — c&X 



(3.27) 

m=iv„ 



This expression evaluated with T from eq. (3.26) is shown in fig. 13 for the four cases in 



- 19 - 



table 3.21. Alternatively, the residue A can be determined by the value of the extrema at 



T 4 (1 + <5 2 ) 4 - 



(3.28) 

t»=t»m,± 

Our first observation is that the results for all three of these expressions seem to agree quite 
well. One might also note how similar the results are for cases V and VII, which both have 
^iim = 0.651. The most significant effect apparent in fig. 13 is that A decreases quite rapidly 
as q increases. This effect is most dramatic in case IV where the results for A are fit well by 
an exponential: A ~ 140 exp[— 0.077 q]. This rapid fall in A with increasing q was the third 
effect which we identified as contributing to the disappearance of the quasiparticle peaks in 
the spectral function. 
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Figure 13: The residue A calculated for the first peak in the four cases given in table 3.21. The 
circles correspond to results for A\ in eq. (3.27), while the squares and diamonds correspond to A^^± 
in (3.28), respectively. 



4. Quasinormal modes 

The basic features of the spectral functions are controlled by the analytic structure of the 
corresponding retarded correlators in the complex frequency plane. Holographically, these 
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poles correspond to quasinormal modes of excitations in the black hole geometry induced on 
the D7-brane [22, 24]. Investigating the quasinormal spectrum for the black hole embeddings 
of interest is a technically challenging problem which we intend to pursue elsewhere [25]. 
However, as described in [9], some qualitative information can be inferred from casting the 
relevant radial equations of motion in the form of a Schrodinger equation. In particular, 
one naively expects that the appearance of quasiparticle peaks with r <C ^ in the spectral 
function would correspond to metastable states supported by a local minimum in the effective 
potential of the Schrodinger equation. Conversely, the absence of a local minima could be 
expected to yield the non-existence of quasiparticles. In this section, we apply this kind of 
qualitative analysis to the transverse vector and pseudoscalar fluctuations of the D7-brane. 
These calculations allow us to gain some qualitative insight into the quasinormal modes for 
these fields - see, e.g., [10, 26]. One might be tempted to apply a WKB analysis in this 
framework to obtain some quantitative results as well, however, as we describe in appendix 
A, such an approach will not yield reliable results. 

Let us begin with (3.11), the equation of motion for the transverse vector, and follow the 
approach described in appendix of D of [9]. Defining 10 Hq = (f 2 FA)/(f 2 G) and Et = hip, 
we find that choosing h = Hq^/F 1 / 2 recasts this equation into 

-d 2 R J + V^ = W 2 ^, (4.1) 



where R* = J°° dp/ y 7 Ho(p) and the effective potential V is given by 



f 

V = —/S. 

P 



-^(F^+q 2 



= Vb + q Vi . (4.2) 

Note that with the above definition for i?*, we have R* — > oo at the horizon (i.e., p — ► 1) 
and R* — ► asymptotically (i.e., p — > oo). We also note that the simple dependence of the 
Schrodinger problem on the frequency and wave-number of the mesons. That is, to 2 plays the 
role of the effective energy while q 2 appears as the coefficient of a new term V\ in the effective 
potential. Further we observe that in the limit that p — > oo, V\ — ► 1 and hence in this UV 
regime, these two contributions can be recombined in the Lorentz invariant form to 2 — q 2 . 

We begin by illustrating the behaviour of Vb in fig. 14. We might compare this to the 
corresponding effective potential discussed for the Minkowski embeddings in [18]. Both there 
and here, the effective potential rises to infinity for large p which simply reflects the infinite 
gravitational potential of the AdS geometry. In the Minkowski case, there is also an infinite 
barrier at p = Rq where the D7-brane embedding ends and so this potential leads to a discrete 
spectrum of (stable) bound states. As described in section 2.1, the D7-branes of interest here 
with d have a narrow neck that extends down to the black hole horizon at p = 1. As a result, 
the infinite barrier at p = Rq above, is reduced to a finite barrier, as illustrated in fig. 14. 



10 Recall that F and G are defined in (3.12) and A, in (2.14). 
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As a result, one's intuition should be that the low lying bound states for the Minkowski 
problem remain essentially unchanged but they are now only metastable because they will 
slowly tunnel out through the potential barrier to the horizon. Further as shown in fig. 14, as 
d increases (with m fixed) , the height of the potential barrier shrinks and so the decay rate 
of these metastable states will increase. Of course, for highly excited states with t» 2 above 
the top of the barrier, the spectrum will be completely changed and in particular, we do not 
expect to speak in terms of metastable states. 




Figure 14: The plot on the left is for Vq vs p for u lim = 0.995 and d = 0.15 (case IV; green), 0.25 
(case I; blue), 0.5 (red). The potential barrier diminishes with increasing d. The plot on the right 
shows the total potential V with q = 10 for the same range of parameters. 

When the meson states also have a finite momentum, we must also account for the 
contribution q 2 V\ to the effective potential. Given the expression in (4.2), one finds that this 
term raises the effective potential by a finite amount are larger values of p. As illustrated on 
the right in fig. 14, while both the local maximum and minimum in the potential are raised, 
the primary effect of increasing q is to raise the minimum relative to the maximum, i.e., the 
barrier shielding the potential well from the horizon is reduced. This effect is also illustrated 
in fig. 15 where we have plotted V — q 2 with changing q - subtracting q 2 ensures that the 
asymptotic form of the resulting curves is fixed, as can be inferred from the discussion below 
(4.2). 

To summarize, recasting the equation of motion in the Schrodinger form (4.1) gives 
an intuitive picture which yields the following key observations: For small q, the effective 
potential has a large barrier, which is responsible for the metastability of the quasiparticles. 
However, this barrier shrinks with increasing q which then explains why the decay rates of 
the metastable states should increase as q increases. Further this barrier actually disappears 
above some momentum which again points to a critical momentum q crit beyond which the 
quasiparticles disappear. 

We can use this analysis to establish a quantitative criterion for this maximum momentum 
q crit . Since our intuition is that metastable states exist because of the potential barrier 
appearing in V, we propose to estimate q crit as the value of the momentum at which the 
barrier disappears and is reduced to a point of inflexion. We have shown these estimates in 
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V-q 2 




Figure 15: V — q 2 vs p for case I (i.e., d = 0.25 and xo = 0.99). The curves correspond to q = 
(blue), 10 (pink), 20 (yellow) and 50 (green). As q increases the barrier decreases and eventually 
disappears. 



fig. 9 for the three cases described in table (3.21). It seems that for cases V, VI and VII 
(with the smaller values of % m ), these estimates are well within the region where the peaks 
are losing their canonical shape and in fact, are quite close to the point where the zero at t»o 
is lifted. The latter was also suggested as a measure of q crit in the previous section. For case 
IV with Vy lm = .995, the result here gives a smaller value of q outside of the region where peak 
is losing its shape and so present calculation seems to underestimate q crit . However, it still 
seems to coincide with the onset of the rapid rise in T. It seems somewhat surprising though 
that the widths could still be relatively small after the barrier in the effective Schrodinger 
potential had disappeared. 

In any event, the advantage of the present approach to estimating q crit is that it is com- 
putationally simple and relatively inexpensive. For example, in fig. 16, we used this approach 
to determine the behaviour of q crit as a function of T/M. As T/M essentially fixes v lim , these 
results are again displaying the strong correlation between q crit and the limiting speed dis- 
cussed in the previous section. In fig. 16, the continuous curves were constructed by fitting 
the data points with the form exp[— 8x]/x 2 . Hence we find in the zero temperature limit, 
q crit oc (M/T) 2 — > co. Furthermore, in the limit that d — > (and hence Tq — > 0), q crit — > oo. 

We close this section with a few technical comments. One may have imagined that the 
decay widths would be calculable in the corresponding Schrodinger problem using a WKB 
approximation. Such a calculation is outlined in appendix A. However, as also explained 
there, this approach does not generally yield reliable results. The essential point is that for 
the WKB approximation to be valid, the change of the momentum over a wavelength must 
be small compared to the momentum itself. Unfortunately, in the examples considered here, 
this condition in not met within the barrier making the WKB analysis unreliable. 
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Figure 16: This plot shows q crit vs T/M for Ex- The blue, red and green points and curves correspond 
to d = 0.25, 0.0005 and 0.0001, respectively. Further the boxes surround the points corresponding to 
cases I, II and III in table (3.19). The continuous curves come from fitting the data points with the 
form exp[— 8x]/x 2 . 



Finally we comment briefly on the effect on the widths from adding angular momentum 
on the internal S 3 . That is, all of mesons considered up to this point have been singlets 
under the internal SU(2) x SU (2) global symmetry of the gauge theory - this symmetry is 
dual to rotations on the D7-brane's internal S s . However, these states only correspond to 
the lowest dimension operators in an infinite family of vector operators transforming in the 
{1/2, £/2) representation of the internal symmetry [18]. So we wish now to consider states 
with nonvanishing I. In fact, we will turn to the pseudoscalar mesons at this point since the 
analysis in this case is somewhat simpler. For non-zero I and q, A p cannot generally be set to 
zero [9] which complicates the analysis for the vector. However, adding angular momentum 
to the pseudoscalar case is straightforward and so we consider this case here. The equation 
of motion for the pseudoscalar s, i.e., fluctuations in cf>, can be written as 



d p (Fd p V) + G(vo 
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q 2 )V + £(£ + 2)HV = 0, 



where the functions F, G, H are defined as 
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It can be shown that the effective Schrodinger potential in this case can be written as [9] 

V = V + q 2 Vx + i(l + 2) V 2 (4.6) 
= -^d p {Fd p h) +q 2 H 2 + £(£ + 2)H 3 , (4.7) 

where h = Hq^/F 1 ^ 2 . Some examples of the effective potential are shown in fig. 17. The key 
feature which one observes from this plot is that introducing £ ^ increases the potential 
barrier. This is not surprising since the D7-brane embedding has a narrow neck near the 
horizon, i.e., the S* 3 is small, and so this increase reflects an angular momentum barrier for 
these modes. We also note that the plots of the potential for the vector meson Et for £ = 
are almost indistinguishable from those for the pseudoscalar even for non-zero q. In fact the 
two potentials differ by less than 1% although their precise functional forms are not identical. 




Figure 17: This plot is for the effective Schrodinger potential for the pseudoscalar field with £ = 0, 1 
and q = 0, d = 0.25, xo = 0.99. The higher potential is for 1 = 1. 



5. Beyond q crit 

As we saw in section 3, the analysis of the spectral functions can only give limited information 
about the dispersion relations. In particular, when the quasiparticles become too unstable 
with increasing q, they no longer contribute characteristic features to the spectral function 
which would allow us to infer f2(q) and T(q). Of course, this simply indicates that the 
corresponding pole in the thermal correlator (3.2) has moved down in the complex plane, 
too far from the real axis to strongly influence the spectral function. The behaviour of these 
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poles can be followed to higher momenta in the present holographic framework with a direct 
analysis of the corresponding quasinormal modes [21, 27, 28]. While we intend to pursue 
this problem elsewhere [25], here we will attempt to provide some qualitative insight into the 
behaviour of the poles at higher momentum and also for the higher quasinormal modes, using 
the intuitive picture of the effective Schrddinger problem given in the previous section. 

Our intuitive picture will also be guided by a comparison with the results for the R 
currents in N = 4 SYM at finite temperature. That is, we start by considering the behaviour 
of the transverse vector modes of a Maxwell field in AdS$. A closely related investigation of 
the longitudinal vector modes and a massless minimally coupled scalar appeared in [27] and 
[28], respectively. 11 The relevant equation of motion for the transverse vector modes Et in 
AdS$ can be found in [9, 21] and can be written as 

„3, Q ^ , 8/ 



d p {p 3 fd p E T ) + - 4 \^, 2 ^ -q 2 jE T = 0, (5.1) 

where / and / are defined above in (2.2). This equation of motion is equivalent to that of 
the analogous modes of the D7-brane gauge field (3.11) with \ = and d = 0. Hence using 
the results of section 4, eq. (5.1) is easily recast into the Schrddinger form 

-clY> + y( P )^ = H>V, (5.2) 



where R* = f^dji/^H^) with H (p) = p*f 2 /8f. The wave-function has been defined as 
V> = Et/H with h = (8p 2 /) -1 / 4 and the effective potential V(p) is given by 



V(p) = Vo +q 2 V 1 = L- 



(5.3) 



/ 2 L32 ^ 8 p 2f 

Again the effective potential in eq. (4.2) reduces to that above upon setting both \ an d 
d to zero. Fig. 18 presents various plots of the effective potential versus the Schrddinger 
coordinate i?* for various values of q. In fig. 20, we also compare this potential (5.3) for the 
supergravity vector modes with that for the analogous D7-brane modes for different values 
of the quark mass and d. 

As noted above, [27, 28] presented an analysis of quasinormal modes for similar super- 
gravity fields in AdS§. Here in our preliminary analysis of the transverse vector, we use 
only the spectral function techniques described in section 3 for the lowest lying modes. The 
analogue of eq. (3.16) is given here by 

P s f P f \f w 2 / 

where S 1S defined in (3.15) with F = p 3 f. The boundary condition for regularity at the 
horizon becomes 5(1) = — 1§ and the spectral function is given by the asymptotic limit in 
(3-17). 

11 There is an erroneous claim in [27] that the results for the transverse vector modes should be precisely the 
same as for that for the massless minimally coupled scalar field given in [28] . 
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Figure 18: The plot shows V — q 2 vs i?*. Blue, red, brown and green are for q = 0,10,20,30 
respectively. The dashed lines indicate the lowest eigenfrequency obtained from the spectral function 
at these q's. 
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Figure 19: The energy f2(q) and width T(q) for the lowest lying mode for the supergravity gauge 
field. The black dashed line on the right corresponds to = q. 



Next we can repeat the analysis of section 3.1 to extract £l(q) and T(q) from the spectral 
function - in particular, our results were derived using the structure of c^X. The behaviour 
of the lowest lying mode are plotted in fig. 19. The most striking features revealed there 
are a) the asymptotic velocity is precisely one and b) the width decreases with increasing q. 
Hence there is a remarkable contrast between this behaviour for the supergravity vector and 
that described for the D7-brane vector in section 3.1. 

We should note that at q = our approach here yields (O(0), T(0)) = (1.51, .85) while 
the precise analytic results are known to be (Q n (0),T n (Q)) = (n,n) [9]. Hence we have 
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an indication of the potential errors in our spectral function results. We should note that 
an analytic expression for spectral function at q = is also known [9] and our numerical 
calculation of spectral function agrees well with this result. However, it is not surprising that 
our technique for extracting the energy and width from 'H is giving imprecise results here since 
we have f2 n (0) = r n (0) and further Af2 ~ T n for the first few poles. We remain confident 
that the qualitative behaviour appearing in fig. 19 is correct. In particular, precisely the 
same behaviour was first found for the supergravity modes studied in [27, 28], e.g., see figures 
3 and 4 in [28]. 

We would like to understand these qualitative features from the structure of the effective 
potential in the Schrodinger equation (5.2). As shown in fig. 18, one of the interesting 
features of V is that it develops a relatively flat plateau in the regime < i?* < 1, which 
corresponds to large values of p. Further, calculating the real part of the effective energy 
Qs = Re(t» 2 ) = Q 2 — T 2 - see appendix A - for the lowest lying level, we find that this 
energy is just above the potential energy at the plateau. If we consider a WKB calculation 
of the corresponding wave-function, our intuition is confirmed about both of the qualitative 
features above. In particular, in this approximation, the amplitude of the wave-function 
would be given by [30] 

H< ' a m> = v^m (5 ' 5) 

where P(R*) is the classical momentum of the particle (with energy Q$) evaluated at the 
position i?*. Hence the amplitude of the wave- function is largest on the plateau and so it 
seems natural that the wave- function has its largest support in this region < i?* < 1. For 
large p, we have -R* ~ 2\/2/ p and so the latter range corresponds to large p. Hence the 
radial profile of the corresponding excitation has its support primarily in the region where 
the geometry is very close to that of AdSs and the redshift effects are minimal. Hence it is 
not surprising that the asymptotic velocity of these excitations is one, i.e., the speed of light. 

Quantitatively we may also observe from fig. 18 that V p i a teau ~ q 2 and so with a small 
width, we naturally have Qs — ^ 2 — 9 2 - Finally while V p i atea u and are rapidly rising 
with q, the value of the potential remains fixed at zero at the horizon which corresponds to 
i?* — > oo. Hence eq. (5.5) indicates that the relative amplitude of the wave-function at the 
horizon is falling \ip\ 2 ~ 1/q. Hence it seems natural that the flux absorbed by the horizon 
should also be decreasing and hence the width of the state would be decreasing. 

In principle, one could verify this intuitive picture with more detailed calculations to 
produce a better understanding of the precise normalization of the wave-function. 12 However, 
we now wish to apply our newly gained intuition to infer the behaviour of the meson dispersion 
relations at very large momentum. In fig. 20, we compare the supergravity potential (5.3) to 
that for the analogous vector modes on the D7-brane for two of the cases (I and V) studied in 
section 3. In both cases, we are in a high-momentum regime where q 3> q crit and so there is no 

12 In the case of near-critical embeddings, a fiat plateau in the potential produced similar effects in the 
quasinormal spectrum at q = [29] . 
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Figure 20: The solid lines correspond to the effective potential in case 1 (left) and V (right) for q = 
80 (blue), 90 (red), 100 (green). The dashed lines correspond to the analogous supergravity potential 
in (5.3). 

evidence of the structure (i.e., a potential barrier and well in front of the horizon) discussed in 
section 4. While the structure of the potentials in fig. 20 differs for > 1 (i.e., deep in the 
bulk of the AdS§ black hole), they are essentially identical for R* < 1 (i.e., in the asymptotic 
region). In particular, the effective potential for the mesons also develops a plateau with 
Vpiateau ~ ^ m this region. Hence it is natural to assume that at very high momentum, the 
poles characterising the thermal correlator for the meson operators will display a behaviour 
very similar to that found in the supergravity analysis. That is, S7(q) should exhibit an 
asymptotic velocity of one at very large q and the widths T(q) should decrease as q increases 
to very large values. 

Hence we are led to conjecture that the full dispersion relations found by studying the 
quasinormal modes dual to the meson operators might take a form as illustrated in fig. 21. 
In particular, the behaviour of the low-lying modes found for q < q crit in section 3 should 
not be indicative of the overall structure. If we consider higher modes, i.e., modes which 
do not lie in the potential well of section 4, we expect that these modes will not exhibit a 
linear dispersion relation with slope v Vlm in the q < q crit regime as the low-lying modes do. 
Rather their dispersion relations will directly approach a linear behaviour in the q S> q crit 
regime with an asymptotic velocity of the speed of light. Similarly their widths are likely to 
be monotonically decreasing as found for the supergravity modes in [27, 28]. We have found 
some evidence of this qualitative behaviour for the higher modes from the spectral functions 
but we do not believe our results are quantitatively precise. There remains the interesting 
question of what should be the high-momentum behaviour of the poles corresponding to the 
low-lying modes, i.e., those which exhibit u lim < 1 for q < q crit . The behaviour suggested 
in fig. 21 is that beyond q crit , S7(q) should also reach t> aS ym = 1 for q > q crit . However, we 
can only speculate on precisely how this should arise. One possibility illustrated in fig. 21 is 
that O(q) should rise up and asymptote to Q = q, which would be similar to the behaviour 
exhibited by the sound mode in [21] - see their fig. 5. This would require passing through a 
regime where dqft > 1, i.e., the group velocity appears to be superluminal, but of course, this 
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need not correspond to a violation of causality - for example, see [31]. Another alternative 
also shown in fig. 21 is that these modes approach d^Q\n m = 1 asymptotically without passing 
through a regime with dq£l > 1. 



n 




Figure 21: This plot shows possible behaviours of the dispersion relations beyond q crit . The lowest 
lying modes (red) exhibit Vn m < 1 below q crit but approach dqQ = 1 in the regime q > q crit . The latter 
limiting behaviour may (dashed curves) or may not (dotted curves) exhibit 'superluminaP speeds, i.e., 
dil/dq > 1. The higher modes (blue) directly approach the limit <9qf2 = 1 in the regime q > q crit . 

We add one further observation on the spectrum of the quasinormal eigenfrequencies. The 
poles corresponding to the quasiparticles discussed in section 4 corresponds to modes with 
support primarily behind the barrier in the effective potential (4.2). Further the discussion 
here also refers to modes with support primarily in the asymptotic region of the D7-brane. 
That is, the support of all of these modes is mainly above the narrow throat in D7-brane 
geometry. However, at least in the regime q < q crit , there should be additional poles cor- 
responding to quasinormal modes whose support is predominantly in front of the potential 
barrier, i.e., in the narrow throat. These would be interpreted excitations of the effective 
string gas that is modeled by the throat [14]. These would be expected to have T ~ f2 and so 
would not lend themselves to a quasiparticle interpretation [4]. It would also be interesting 
to understand the behaviour of these quasinormal modes and their interplay with the modes 
with support above the throat as we move into the regime q > q crit . Certainly, as shown in 
fig. 20, the structure of the effective potential differs from that in the supergravity problem 
for large R* (small p), hence it could be that it still supports an additional set of modes in 
this regime with distinct physical characteristics. 
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6. Discussion 

In this paper, we have used holographic techniques to calculate the spectral functions and to 
study the dispersion relations of meson quasiparticles moving through a thermal plasma of 
a strongly coupled N = 2 super- Yang-Mills theory. The quasiparticle peaks in the spectral 
functions arise from poles in the corresponding thermal correlator (3.2) at 0, — iT in the lower 
half of the complex frequency plane. Considering the energy f2 and the width V as functions 
of the spatial momentum q, the broad picture which emerged was that there were two distinct 
regimes distinguished by a critical momentum q crit . In the low-momentum regime q < q crit , the 
spectral function exhibited clear quasiparticle peaks. By following their position and shape 
with growing q, we were able to estimate the the dispersion relations, f2(q) and T(q), for 
the low- lying resonances. In the high-momentum regime q > q crit , the poles have moved too 
far into the complex plane and no quasiparticle peaks are evident in the spectral functions. 
However, we examined the effective Schrodinger equation governing the dual quasinormal 
modes and provided some qualitative insight into the dispersion relations in this regime. 

One of our key results for the q < q crit regime is that the quasiparticles approach the 
same limiting velocity found for the case of stable mesons [7, 12]. This result emphasizes 
that this limiting velocity will be a universal feature in any holographic model. While the 
precise form of the limit may depend on the particular model, it arises from gravitational 
redshifting in the background geometry, as indicated by (1.1). Hence it will apply for any 
gauge theory excitations that have a dual description in terms of radially localized modes in 
the dual geometry. So, for example, this robust feature would also apply in the adjoint sector 
to 'glueball' excitations that are dual to a wavepacket of supergravity modes localized in the 
radial direction. 

Of course, not all excitations of interest in a holographic model need not be radially 
localized. For example, massive quarks in the N = 2 super- Yang-Mills theory (with n q = 0) 
are represented by extended strings stretching down from the D7-brane to the horizon in the 
dual gravitational description. Of course, the same gravitational redshift is observed to have 
interesting physical effects when these quarks are in motion - e.g., see [32]. Further, using 
holographic Wilson lines, a similar effect was observed to lead to the dissociation of a heavy 
quark bound state moving at a finite velocity through a strongly coupled N = 4 plasma [33] . 

The second interesting observation for the low-lying mesons in the q < q crit regime is that 
their widths show a dramatic increase as the momentum approaches q crit . In our holographic 
gravity model, the rise in V is easy to understand. As observed for the Minkowski embeddings 
[7], with finite q, the excitations on the probe D7-brane feel an extra potential which pushes the 
support of the radial profile of the mesons down to minimum radius on the brane. Intuitively, 
this can be understood as when the momentum is introduced, the brane fields carry extra local 
energy density and so feel a stronger gravitational force pulling them towards the horizon. 
As described in section 2.1, we have destabilized the mesons by introducing a narrow throat 
which extends down to the horizon from the point where the brane would have otherwise 
closed off. Hence the additional potential due to the momentum naturally pushes the meson 
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fields down this throat making these states decay more quickly. This intuitive picture can be 
made more quantitative with the Schrodinger framework introduced in section 4. The latter 
approach also made clear that the potential barrier, which maintained the metastability of 
the quasiparticles, vanishes above some momentum. The disappearance of the barrier in the 
effective potential provides a precise mathematical criterion with which to define q crit . While 
this is a natural definition, we should add that the distinguishing momentum q crit remains a 
rather qualitative concept and another definition was also considered in section 3. 

In the high-momentum regime, our qualitative investigation of the effective Schrodinger 
equation suggests that the asymptotic velocity is v asym = 1 for q » q crit . Here the modes 
should be radially localized but evade having a limiting velocity less than the speed of light 
because they are not localized near the black hole horizon. Rather, our suggestion is that 
as q increases to very large values, the support of quasinormal modes becomes increasingly 
focussed at very large radius, i.e., towards the AdS boundary, where the gravitational redshift 
becomes vanishingly small. This conjecture is largely based on making an analogy with similar 
supergravity fields for which some explicit results for the quasinormal frequencies are known 
[27, 28]. The latter results further suggest that in this high momentum regime the widths 
should decrease with increasing momentum. To clarify the details of the behaviour in this 
high-momentum regime, one would have to examine the quasinormal modes directly [25]. It 
would also be interesting to better understand the physical differences between the present 
case and that with n q = in which the mesons are display u lim < 1 for arbitrarily large q. 

Still given the present understanding, an interesting physical picture has emerged. Namely, 
in the low-momentum regime, the quasiparticles are strongly coupled to the deconfined plasma 
of the adjoint fields. The gravitational redshift leading to v lim < 1 is the geometric description 
of this strong coupling. An interesting problem would be to map the 'glue' cloud associated 
with these moving mesons, e.g., along the lines of [34]. We expect that the energy density of 
this halo must be rapidly increasing to maintain the meson's velocity at v Um < 1. However, in 
the high-momentum regime, the quasiparticles and the adjoint plasma are no longer strongly 
coupled to each other and the meson excitations can achieve the speed of light. Of course, 
the behaviour in this high-momentum regime restores the intuitive picture that one might 
acquire from considering the theory at weak coupling, namely, that at high momentum, the 
quasiparticles should be largely unaffected by the surrounding plasma. 

In the low-momentum regime, the methods which we applied in section 3 allowed us to 
estimate not only the positions of the poles (from and T) but also their residues A. The 
latter were found to decrease with increasing q, as illustrated in fig. 13. As noted in [35], the 
parameter dependence of A can be important in determining the overall form of the spectral 
function and clearly here the fall in the residues plays a role in the disappearance of the 
resonances in the spectral functions. However, one might note these residues are indicative of 
the coupling of the relevant operator to the thermal bath, i.e., of how effective the operator is 
in generating the relevant quasiparticle excitations. Hence for the present purposes, we need 
not think of A as a physical characteristic of the quasiparticle itself. 13 In contrast, f2(q) and 
13 However, this coupling may become physically relevant for a certain quasiparticle decay channel if in the 
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r(q) certainly characterize the basic physical behaviour of the quasiparticles. In particular, 
in order for a certain pole in the thermal correlator to be considered a quasiparticle in the 
first place, it must satisfy the Landau criterion r <C fl. Intuitively, we may understand this 
requirement by observing that the corresponding wave function has a factor exp[— iClt — Tt]. 
Hence with r <C £1, there are a large number of oscillations before the excitation is damped 
out. 

In section 3, we found another important reason that the quasiparticle peaks were quickly 
absorbed into the background is that the spacing between neighbouring poles shrinks with 
growing momentum, as is illustrated by the dispersion curves in fig. 6. Certainly at a practical 
level identifying the peaks in the spectral functions requires T <C A$7 where AQ is the spacing 
between neighbouring poles (along the real axis). As illustrated with fig. 10 for case VI, the 
peaks are begin to coalesce at roughly the point where they vanish from the spectral function. 
At this point, the poles of interest are not isolated and rather there are a large number of 
poles with roughly equal spacings and similarly growing widths, which produces some smooth 
continuum (rather than individual peaks). A quasiparticle interpretation is inappropriate 
then. Instead the infinite collection of poles is collectively generating the smooth background 
which remains in the spectral function. 

It is interesting to consider how well the peaks of the spectral function were reconstructed 
given the estimates we produced for the relevant parameters in section 3. Fig. 22 shows a 
typical peak which corresponds to case IV with q = 50. Also shown is the Breit-Wigner peak 
(3.3) reconstructed with various estimates of Q, T and A, produced by examining <9j 'H as 
described in section 3. While this reconstruction matches the shape of the peak near the 
maximum quite well, there are two other notable features which are quite apparent. First the 
reconstructed peak appears displaced slightly towards smaller to relative to where = 0. 
Of course, this is to be expected since there is a rising background that also contributes to 
the spectral function. The second feature, which we found surprising, is that this background 
seems to be quite small. That is, the maximum of reconstructed peak matches very well 
with the maximum value of the spectral function. More generally, we found that using 
our parameter estimates to reconstruct the quasiparticle peaks gave a background that was 
typically only between 10 and 20% under the top of the peak. 

Given the limitations in studying the spectral functions, it would also be very interesting 
to investigate the quasinormal modes directly, following e.g., [10, 26]. While this presents 
additional technical challenges, it would allow for a detailed verification of the quasiparti- 
cle dispersion relations found here in the low-momentum regime but also to establish the 
behaviour conjectured for the high-momentum regime [25]. At this stage, we would like to 
make the following observation. It is possible to extract the quasinormal modes directly from 
the spectral function. If to = — iT is a quasinormal mode, this will show up as a pole in the 
spectral function Ti. Since we have an approximate idea about the location of the peaks from 
the spectral functions as well as an idea about the widths, we have a reasonable initial guess 

underlying theory the operator couples to other physical fields, e.g., photons [36]. 
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50.0 50.1 50.2 50.3 50.4 

Figure 22: The red curve shows the spectral function for case IV with q = 50. The blue region shows 
the peak (3.22) as reconstructed with the parameters Q, T and A. calculated in section 3. The center 
(dashed) line in the shaded region uses (fl,T,A) — (t»o,r, A\). The bottom (solid) boundary of the 
shaded region is generated with (Cl,T, A) = (tt>o,f + AT,A 2 ,+), while the top (thick solid) boundary 
uses (il, T, A) = (to , f - Ar, A 2 _). 

which can be fine-tuned to locate the pole. It turns out that for larger q's, the fine-tuning 
required increases. The following table is a comparison between this method and the method 
of estimating parameters from for case VI. 
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As is clear from the table, for larger q the quasinormal mode analysis predicts a lower V. 

While we see a dramatic rise in the width T of quasiparticles in the present model, it 
remains to understand if this increase is a universal feature which emerges in any holographic 
model. Of course, an even more important question, is whether or not such an effect is realized 
in the strongly coupled quark-gluon plasma of QCD. Certainly both of these questions deserve 
further investigation. 

At present, investigating spectral functions with non-zero momentum is an important 
direction of ongoing research [37, 38]. Only partial results are available since the existing 
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methods seem to be inadequate for this problem. Recent results, which generally use the so- 
called potential approach, indicate the spectral functions flatten with increasing momentum 
[37, 38] - an effect consistent with our present findings. 

Of course, if these effects found here in our holographic studies are realized in QCD, 
they would have interesting implications for experiments at RHIC and LHC. In particular, 
as suggested in [12, 33], they could lead to a significant additional suppression of J/* or 
other heavy quark mesons with large transverse momenta. Certainly a critical momentum 
where the quasiparticle widths rise dramatically would produce a dramatic effect. Another 
potentially interesting effect was outlined in [36]. One distinct feature resulting from a limiting 
velocity t> lim < 1 is that the four momenta of quasiparticles become null at some point. For 
example, if we drop the higher order terms in the asymptotic dispersion relation (3.20), then 
the quasiparticle's dispersion relation crosses the null cone at 



Following the holographic techniques of [39] , one then finds that this crossing of the null cone 
produces a peak in photon production from charged quasiparticles [36]. Such a resonance 
would then be a distinctive experimental signature of v Um < 1. Of course, as speculated in 
fig. 21, the dispersion relations may cross the null cone twice and exhibit a regime dqtl > 1, 
which may produce further dramatic signals. 

From (6.1), we observe that this enhancement is pushed to infinite momentum as t> lim 
approaches one. Hence, if v liui is still close to one at RHIC or LHC, this signal would only 
appear at very large momenta. Further, if the quasiparticle width increases too quickly, this 
enhancement in the photon production would likely be washed out. The point where the 
quasiparticle 'four-momentum' is null is indicated in fig. 9 for the four cases in table (3.21). 
For cases IV, V and VI, one indeed finds that this null momentum seems to lie in the regime 
where the width is growing rapidly, while case VII seems to present an exception to this 
rule. Further as can be seen in fig. 7 for the more stable cases presented in table (3.19), the 
quasiparticle dispersion relations seem to cross the null cone well away from the maximum 
momentum. Hence it seems this question also requires further study. 
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A. WKB approximation 



In this section, we outline how WKB calculations might be applied to produce an approximate 
value for the quasinormal eigenfrequencies. Recall that these eigenfrequencies are determined 
by solving the relevant wave equation with boundary conditions of an ingoing wave at the 
horizon and of only the normalizable mode asymptotically [10]. In (4.1), the relevant wave 
equation for the transverse vector modes has been caste into the form of a one-dimensional 
Schrodinger with an effective energy 



E e ff = fis - i^s 



(A.l) 



= t» 2 = n 2 - r 2 - 2i nr . 



The quasinormal eigenfrequencies correspond to the positions of the poles in the thermal 
correlators and so we denote the eigenfrequency to = Q — iT above in keeping with the 
notation of section 3. 

In a typical case of interest, the effective potential has a form as illustrated in figs. 14 
or 17, with a well that is separated from the horizon at p = 1 by a large potential barrier. 
The relevant modes corresponding to metastable meson states are then bound states with 
support primarily in the potential well but which slowly tunnel out through the barrier. The 
idea then is to approximately determine E e ff for these modes using WKB techniques - see, 
e.g., [30]. 

First the real part of the effective energy Q$ is obtained by fine-tuning Qs so that 

rRa 

(A.2) 
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where #2,-^3 are the classical turning points in the potential well where V = fls- The first 
peak in one of the spectral functions would correspond to n = 1. Next within the WKB 
approximation, the decay rate T$ is defined as [30] 

-R 2 
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Here R±, R2 denote classical turning points in the barrier and implicitly fig denotes the WKB 
energy eigenvalue determined by (A.2). The normalization constant N is given by 
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Given these WKB results for the Schrodinger problem, we may use (A.l) to determine 
the eigenfrequency for the corresponding quasinormal mode: 
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Figure 23: The validity of the WKB approximation can be tested by considering k'/k 2 versus p. 
These are plotted here inside the barrier for d = 0.25, xo = 0.99. At the turning points a deviation 
from unity is expected but significant deviations are seen well away from those points, making the 
WKB approximation (A. 3) unreliable. 



In the second of each of these expressions, we have presented the leading terms in an expansion 
with Tg/fis <C 1, which is a necessary condition for the validity of the WKB approximation. 
The latter calculation focussed on the well and barrier appearing in the effective potential 
V(R*) while completely ignoring the horizon. In the results of, e.g., [9, 10, 21], one sees that 
for an effective potential is a monotonically rising function T ~ or Tg > fig. Clearly, the 
horizon has an important effect in determining the eigenfrequencies and, in particular, the 
decay width in this situation. Hence one cannot expect that the WKB calculations outlined 
above will be reliable in this regime. In particular, when the barrier is getting small, the 
WKB calculation would yield an tunnelling rate that becomes large but the precise value of 
Ts would depend crucially on the structure of the potential near and the boundary conditions 
at the horizon. 

However, as discussed in section 4, there are additional subtleties in the WKB approx- 
imation. For the WKB calculation of the tunnelling rate to be valid, the change of the 
momentum over a wavelength must be small compared to the momentum itself. That is, if 
we define k = y/V — to 2 under the barrier, then we must require [30], 

(A.6) 



k(R*) 



Unfortunately, in the examples considered in the paper, this condition is typically not 
satisfied for a large region within the barrier, essentially because the shoulders of the barrier 
are too steep - e.g., see fig. 23. This makes the WKB calculation (A. 3) of Tg unreliable, 
in particular in the regime where T begins to increase dramatically. We have done WKB 
calculations of T in various cases when T/Tq ~ 1 and found the WKB results agreed with 
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those in section 3 to within roughly 15% despite the issues discussed above. Note that the 
condition (A. 6) is quite generally satisfied for the calculation (A. 2) of fig, the real part of 
the effective energy. Hence in the regime Ts/^s <S 1> the WKB approach would still provide 

1/2 

an accurate approximation for 17 ~ 0,^ . Further, one can still gain some qualitative insight 
into T by considering when the barrier disappears in the effective potential, as discussed in 
section 4. 
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